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ABSTRACT 

A nonlinear representation for the rapid-pressure correlation appearing in the Reynolds stress equa- 
tions, consistent with the Taylor-Proudman theorem, is presented. The representation insures that 
the modeled second-order equations are frame-invariant with respect to rotation when the flow is 
two-dimensional in planes perpendicular to the axis of rotation. The representation satisfies real- 
izability in a new way: a special ansatz is used to obtain, analytically, the values of coefficients 
valid away from the realizability limit: the model coefficients are functions of the state of the tur- 
bulence that are valid for all states of the mechanical turbulence attaining their constant limiting 
values only when the limit state is achieved. Utilization of all the mathematical constraints are 
not enough to specify all the coefficients in the model. The unspecified coefficients appear as free 
parameters which are used to insure that the representation is asymptotically consistent with the 
known equilibrium states of a homogeneous sheared turbulence. This is done by insuring that the 
modeled evolution equations have the same fixed points as those obtained from computer and labo- 
ratory experiments for the homogeneous shear. Results of computations of the homogeneous shear, 
with and without rotation, and with stabilizing and destabilizing curvature, are shown. Results 
are consistently better, in a wide class of flows which the model not been calibrated, than those 
obtained with other nonlinear models. 
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1. Introduction 

Most turbulence models are devised for use in inertial coordinate systems. Ad hoc changes are 
then made to reflect the unusual effects on turbulence of swirl, curved streamlines or rotation of 
the coordinate system: one makes the length scale depend on Richardson number or adds terms 
to the dissipation equation calibrated to observed behavior. This approach does not make use of 
the mathematical requirements that the dependent variables and their evolution equations must 
satisfy and leads to models that perform poorly when used in situations substantially different 
from the benchmark flows for which they have been calibrated. There is no reason a second-order 
modeling method cannot be successfully applied to a high Reynolds number rotating turbulence: 
no new unknown terms appear in the equations and, for moderate Rossby number, very little of 
the phenomenology on which the technique is based changes. 

The effect of rotation on the second-order moments is felt through the rapid-pressure-velocity cor- 
relation and the Coriolis terms. The difficulty with the equations, as modeled presently, can be seen 
when the equations are transformed to a rotating coordinate system: they are not materially-frame- 
indifferent in the two-dimensional limit, Speziale (1981), Hide (1977). This phenomenon, predicted 
by the Taylor-Proudman theorem requires that when the velocity field is two-dimensionalized by 
rapid rotation, without components along the axis of rotation, the equations must be independent 
of rotation. The problem with the current modeled second-order equations results from the inabil- 
ity of the rapid-pressure correlation model to reflect the physics embodied in the Taylor-Proudman 
theorem. This is another form of realizability: one must obtain, for a bounded flow, in the limit of 
rapid rotation, a frame-indifferent turbulence. In such a turbulence the rapid-pressure correlation 
appearing on the right hand side of the second order moment evolution equations equals the Cori- 
olis terms appearing on the left hand side leaving the equations frame- in different. A rapid-pressure 
model consistent with these facts, reflecting more of the information contained in the Navier-Stokes 
equations, is required. 

As a general tool is being developed to compute a wide class of flows for which there may not be any 
well documented benchmark flows with which to calibrate coefficients, it is necessary to construct a 
model from first principles, incorporating more of the physics. In the present new representation for 
the rapid-pressure strain correlation the use of all constants that do not come from first principles 
have been minimized: This is done this by requiring the representation to have the proper behavior 
in five different limits: 1) frame invariant with respect to rotation when the eigenvalue of the 
Reynolds stress tensor, < U{Uj >, along the axis of rotation vanishes, 2) the realizable limit in 
which an arbitrary eigenvalue of the Reynolds stress vanishes, 3) the joint-realizable limit in which 
an eigenvalue of the tensor < 89 >< U{Uj > — < 9u{ >< &Uj > vanishes, 4) the isotropic limit in 
which the anisotropy tensor, = 0, and 5) asymptotic consistency with a stationary state of the 
turbulence in which D/Dt = 0. Of these five principles the last one is a statement based on 
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experimental observation rather than mathematical fact. It is tacitly asusmed that for a specific 
class of flows, to a reasonably suitable approximation, there exists an equilibrium state to which 
the flow relaxes upon the removal of any disturbing forces. 

Assuring the proper frame-invariance of the modeled second moment equations is done by requiring 
that in the limit of a turbulence two-dimensional in planes perpendicular to the axes of rotation 
the rapid-pressure velocity correlations in the heat-flux and Reynolds stress equations satisfy the 
’’geostrophic” constraints, Ristorcelli (1987) and Ristorcelli and Lumley (1991b): 

^pqn^n^ipqj = fqjn^n ^ ^ 

tpqn^nXpqi ~ €q{ n $l n <C Uq6 > . 

The term geostrophic is borrowed from the meteorological literature where it is used to describe 
the low Rossby number balance between Coriolis and pressure forces in the evolution equations for 
large scale weather. The constraint, however, is independent of how the flow is two-dimensionalized 
and is, therefore, independent of Rossby number. The requirements of the Taylor- Proudman 
theorem are subsumed by the requirements of the principle of two-dimensional frame invariance 
(2DMFI) first put forward, in the context of turbulence modeling, by Speziale (1981). A flow that 
is two-dimensionalized due to stable stratification or magnetic forces is also frame-indifferent in 
the two-dimensional two-component limit and the rapid-pressure correlation must still satisfy the 
geostrophic constraint for arbitrary rotation. 

The requirement of two-dimensional material-frame-invariance comes from first principles: it is 
a rigorous limit of the Navier-Stokes equations for the fluctuating velocity, Speziale (1990). The 
principle of 2DMFI seems to have been given short shrift. Reynolds (1989), in a rapid distortion 
theory analysis of a decaying turbulence, has shown that rotation has a distinct effect on the sec- 
ond invariant, II, of the anisotropy tensor, 6 t j, and goes on further to show that the equations for 
the invariants using the current rapid-pressure models are independent of rotation. He concludes 
that such models, because they don’t show an explicit dependence on the rotation in the evolu- 
tion equations for the invariants, are inadequate for rotating flows. A reference to the evolution 
equations for the invariants, for example Speziale et al. (1991), shows that both the production 
and dissipation, which are functions of the rotation through their dependence on the Reynolds 
stresses, appear in the evolution equations for the invariants. The equations, except for the case of 
a decaying turbulence in which there are no production terms, are not independent of rotation. 

Shih and Shabbir (1990), following Reynolds argument decide that since the rotation terms do not 
explicitly appear in the evolution equation for the invariants (which is true for any rapid-pressure 
model that is a function of symmetric tensors) that the use of a properly frame- in different model will 
not produce the proper turbulence behavior. This argument is specious: the same argument applied 
to the trace of the Reynolds stress, the kinetic energy of the turbulence, q 2 =< u p u p >, one would 
abandon refinement of rapid-pressure models because they do not affect the energy of the turbulence 
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- the evolution equation for q 2 being independent of the rapid-pressure modeling. However, different 
rapid-pressure models produce different Reynolds stresses and different turbulence energies in the 
same flows. Even if the models fails to account for some predictions of RDT for a restricted class 
of flows neglecting results that are rigorous consequences of the Navier-Stokes equations, Speziale 
(1990), are not justified. (It should be kept in mind that the RDT limit violates various equilibrium 
assumptions made in the mathematical development of the second-order modeling method). The 
2DMFI constraint is relevant to flows other than rotating flows: any flows in which the largest 
scales of the motion are two- dimen sionalized by strong body forces will come arbitrarily close to 
being rotationally frame invariant. A dramatic example of this might be the collapse of the wake 
in a stably stratified environment. Shih and Shabbir (1990) also cite Reynolds (1987) as having 
shown that the 2DMFI constraint does not provide any new information. Here it is found that 
the geostrophic constraint produces six additional linearly independent constraint equations for the 
rapid-pressure model. 

The present article is a description of the derivation and validation of a rapid-pressure model 
that satisfies the principle of material-frame-indifference in the two-dimensional limit. The present 
derivation of the 2DMFI model has been given earlier, in departmental reports in Ristorcelli (1987, 
1991) or Ristorcelli and Lumley (1991b). In those previous developments of the 2DMFI model 
several free parameters were, for simplicity, set to zero. In this paper the stationary points of the 
homogeneous shear, as reflected in the structural equilibrium assumption, Dj Dt b{j = 0 where 
bij =< U{Uj > / < u p u p > — 1/3 6ij, are used to set the unknown free parameters occurring in the 
model. The algebraic equations describing the stationary states of the homogeneous shear along 
with the experimentally determined values of the stationary state are used as constraints to set 
the free parameters in the model. This insures that the fixed points of the modeled differential 
equations are the same as the fixed points of a particular physical flow. This calibration is done 
without sacrificing any of the mathematical principles built into the model. 

The next section of this article defines the problem. A third section presents a derivation of the 
tensor polynomial model for the rapid-pressure correlation, A discussion of the constraints as well 
as a derivation of the new geostrophic constraint is given. A new ansatz for the rapid-pressure 
correlation produces a model valid away from both the geostrophic (two-dimensional in planes 
perpendicular to the axis of rotation) and realizability limits (two-dimensional in that one of the 
eigenvalues of the Reynolds stress tensor disappears). The issues of realizability in the light of 
some recent work, Speziale (1993), are discussed. A subsequent section uses results from matrix 
algebra to collapse the model to a structure similar to other rapid-pressure models with some 
interesting differences. The model is seen to have the same tensor bases as the FLT model and in 
the limit of a planar mean flow it has the same tensor bases as the quasi-linear SSG model. The 
two subsequent sections show computations done with the model and compare the results to other 
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models. Following this a general discussion of modeling issues, in the context of realizability and 
the use of stationary states as the only self-consistent asymptotic limit for calibrating turbulence 
models is given. 

2. The rapid-pressure correlation 

In an incompressible turbulence in a rotating coordinate system, with buoyancy effects in the 
Boussinesq approximation, the Reynolds stress equations have the following form: 

D/Dt < UiUj > +2(tikp < UpUj > < u v U{ >)Sl k Ro~ l = + < > (3j+ < 6uj > fa 

[<C UjUp ,p f 'Ui'Up Uj ,p ] UiUjUp ,p 

-[< P,j Ui > + < Uj >] + Re~ l < UiUj >,pp - 2Pe _1 < Uj, p > 

> +2( pt kS l k < 0u p > Ro " l = -[< Ouj > Ui t j+ < U{U 3 > Tj]+ < 99 > /?,■ 

- < > i j — < p t i0 > 

+Re~ l (l + Pr” 1 )(< 0Ui >jj -2 < 0 y3 u ii3 > . 

The velocity has been normalized by a characteristic velocity u c and the Rossby number is Ro = 
u c /QR c where R c is a length scale and f) the rotation rate of the frame of reference. The gravity 
and rotation vectors are aligned with the 3 axis. Our concern is with the pressure-velocity and 
-temperature correlations, < p>jUi > and < p^9 >. An equation for the pressure fluctuations 
comes from the divergence of the Navier-Stokes equations for the fluctuating velocity 

Xl\it H“ XI j U i + U j U| XL^Xlj ^ ,j ~\~2€i kp ^l k U p Ro — “I" R^ 

which produces a Poisson equation for fluctuating pressure. The standard linear decomposition 
recognizes three terms 

-p r t H = 2 [I/,-, p +e p ikttkRo~ 1 ]u p ,i 

~P s a = u > v u j >. ~ < «i « >«■ > 

P,ii = ft*,. 

where p r ,p a ,p b are respectively the rapid-pressure, the slow or return to isotropy pressure, and the 
buoyancy-pressure. The effects of rotation are felt through the rapid-pressure, p r . Solution of the 
Poisson equation for the rapid-pressure is by application of Green’s theorem 

= - (4jr) _1 - x'). 

It is the moments of the solution that are required to close the second-order equations. For a 
homogeneous mean field, more than an integral scale away from any surfaces, a staightforward 
interchange of the order integration and averaging produces: 

^ p T ^ 4" < P n XI j > — —2[Uq ,p -^-€ pqk H k Ro ][A^pgj 4" ^ jpqi] 

< P r n 6 >= -2[f/ 9 ,p +€p qk il k Ro~ l }X piq 

where : 

Xpi q = (4ff) -1 / < 0(x)u p (x') >,;y dx 7(x - x') 

Xtpqj = (4 tt) - 1 / < u,(x)u p (x') > vV dx'/(x - x'). 
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The construction of a tensor polynomial model for the volume integrals of the two-point correlations, 
Xijki and Xijk, is the subject of this article. The rapid-pressure covariance integral appearing in 
the heat flux equations, Xij ky is treated as the Reynolds stresses and the heat fluxes are linked 
through Cauchy- Schwarz inequalities, < 98 >< U{Uj > — < 9ui >< Ouj > > 0 and their modeling 
cannot be done independently. The rapid-pressure term appearing in the Reynolds stress equations 
necessary to treat the mechanical turbulence problem is the focus of this article. 

3* Obtaining a representation for the rapid-pressure correlation 

A constitutive relation in which the integral of the two-point covariance is parameterized by a local 
function of the anisotropy tensor and heat-flux vector is proposed. The most general forms of such 
relationships are the following tensor polynomials 

Xijki/ < Up%> >= AxSijSki + A 2 (8i k bji + fiirfjk) 

+ A 3 6ijbki + A 4 bij6kt + A 5 (bik6ji + bu6j k + $ikbji + 6ubj k ) 

+ AeSijblt + Art^Ski + + b h 6 j * + + M?*) 

+ Agbijbki + A\o(bikbji + bubj k ) 

+ A\\bijb\i + Awbijbki + Ai 3 (tif k bjt + b^bjk + b^b 2 ^ + 6 * ib 2 k ) 

+ Aubtyh + Ais(b? k b? l + b} x b 2 j k ) 

Xpkj — D\ < Ov,p > Sfcj -f* D 2 {< 0u k ^ ^ ^ ^pfc) 

+ D s < 8u v > b kj + D 4 (< 0u k > 6 p j+ < Ouj > b pk ) 

+ A> < 6u p > b 2 kj + De(< 8u k > b^+ < Ouj > b 2 k ) 

+ [D^bqpbkj + D$(bq k 8pj + bqj6p k y\ < 9tlq > 

+ [Z?9 bqpb k j + D\${bq k bpj + bqjbp k y\ <C 0U q > 

+ [Dnb qp b 2 kj + D X2 {b qk b 2 - 4 - b q jb 2 k )] < 9u q > 

+ [Di 3 b 2 6 kj + D l4 (b 2 6 P3 + 6 y p *)] < 0u g > 

+ [Disbqphj + D\e(b k b p j + b jb pk )] < 9u q > 

+ [Dvr b qp b 2 kj + Di S (b qk b 2 vj + b qj bl k )] < 9u q > 

where b{ 3 is the anisotropy tensor bij =< UiUj > / < u q u q > —1/3 6{j and < 8u{ > is the turbulent 
heat-flux. Following Pope’s linearity principle only terms linear in the heat-flux are kept. The A{ 
and Di are functions of the invariants of b^ and < 9ui >. 

Parameterizing the integral of a two-point correlation in terms of the local anisotropy tensor is a 
substantial simplification requiring consideration. For a turbulence with short term memory and 
limited awareness Lumley (1970) has discussed the conditions under which such a constitutive 
relation is tenable and carried out a similar expansion procedure, Lumley (1967), indicating how 
the truncation errors scale. ^From one point of view, the constitutive relation proposed can be seen 
as the first term in a functional Taylor series expansion for the rapid-pressure correlation. As the 
correlation decays with distance the primary contribution to the integral will come from regions 
within an integral length scale of the local position - thus, in a homogeneous turbulence, the first 
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term will constitute a good approximation. Retaining higher order terms of the functional Taylor 
series expansion which involves spatial and temporal derivatives of substantially complicates 
the problem. It is expected that the retention of only the first order term captures enough of the 
physics to allow prediction suitable for engineering purposes. 

In applying the above constitutive relation in non-inertia! systems there are several phenomena, 
peculiar to rotating flows, that may be in conflict with the assumptions of the turbulence closure. 
In particular the inertial wave field associated with the rotation may interfere with: 1) the energy 
cascade from the large to the small scales of the flow, 2) the universal equilibrium assumed for 
the small scales of the flow, and 3) the assumed steadiness of the mean flow. These issues require 
consideration. 

In a mean flow rotating with speed fi there will be inertial oscillations with frequency less than or 
equal to 2f2. For a quasi-steady assumption to be valid, changes in the mean flow must be slow with 
respect to the turbulence’s ability to adjust to the imposed changes. This requires that 2fi < e/k 
and becomes a lower bound for a Rossby number, Ro = s/2Q,k > 1, for which an equilibrium theory 
is appropriate. 

Similar to a stably stratified density field in which stratification inhibits particle motions in the 
vertical direction, rotation inhibits transverse displacements of fluid particles. Similar to the radius 
of gyration of a charged particle in a magnetic field, the lateral displacement of a turbulent fluid 
element in a rotatingframe can be characterized by alength scale (q 2 / 3) 1 ^ 2 /2fl where q 2 =< UjUj >. 
This length scale must be larger than the turbulence length scale £ ~ (g 2 /3) 3 ^ 2 /e to insure that the 
transverse confinement by the Coriolis forces does not affect the fluctuating field. This produces a 
similar bound on the turbulent Rossby number Ro — e/2$lk >3/2. 

The phase coherence necessary for the cascade of energy to the smaller scales of the motion will 
be interfered with if the production scales of the motion k£ ~ 1, where k is the wavenumber, are 
subject to an inertial wave field. A different Rossby number can be defined as a ratio of the vorticity 
of the production scales of the motion to the background vorticity Rot = (q 2 /3) l ^ 2 /2il£ = 3e/2 $lq 2 
using e = (q 2 /3) l ^ 2 /£. A spectral Rossby number can also be defined as Ro(k) = u(k)/2Q,£(k) = 
(kE(k)) 1 / 2 /2Q,(2n / k) which using the inertial range scaling E(k) = ae 2 ^ 3 K~ 5 ^ 3 and e = (q 2 / 3) l l 2 /£ 
becomes Ro( *) ~ 0 .2 (K£) 2 / 3 Ro t . The effects of rotation decrease as the wave number increases. For 
the inertial oscillations associated with the rotation not to interfere with the cascade mechanism 
Ro(k) > 1 for k£ ~ 1 is required. Thus for Rot > 5 the usual parameterization of the spectral 
cascade rate, £, in terms of the energy containing scales of the motion is appropriate. For Ro t < 5 the 
current dissipation equation requires modification. How the dissipation equation is to be changed 
to account for the effects of rotation on the cascade rate is an unresolved issue that is the topic 
of current research. It is, however, clear that the assumption of the small scale equilibrium with 
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the large scales of the motion is valid in most high Reynolds number rotating flows of interest: 
the Rossby number of the dissipation scales of the motion is Ro e = Re/ Rot and therefore the 
dissipation scales are insensitive to the effects of rotation. 

Having dealt with some of the phenomenological issues, the mathematical constraints required of 
the pressure- velocity correlation are now investigated. 

There are several physical inequalities and mathematical identities the X must satisfy. These prin- 
ciples are used to obtain a set of algebraic constraint equations for the Ai and D t . For an arbitrary 
three-dimensional turbulence the tensor polynomials must satisfy the symmetry constraints, 

Xijkl - Xijik 

Xi ik i = Xjiki 

Xijk = Xikj. 

These symmetry constraints are built in to the assumed form of the tensor polynomials. For an 
arbitrary three-dimensional turbulence the tensor polynomials must also satisfy the constraints of 
normalization, continuity: 

Xijkk = < u x Uj > X x kk = < Ui9 > 

Xijjk — 0 X%j{ — 0. 

Note that a contraction of the integral of a two-point statistic is a local one-point statistic. 

The tensors < U{Uj > and < 90 >< U{Uj > - < 0u x >< Ouj > are positive semi-definite. This 
reflects the fact that the energy of the turbulence is always positive and that the magnitude of the 
correlation coefficients between the various components of the tensors be bounded by one. These 
facts lead to the ’’realizability” and ’’joint-realizability” constraints which specify the behavior of 
the correlations when specific limit states are approached. The relevant portion of the Reynolds 
stress transport equations, in principal axes, requires that 

D / Dt ^ rs - / [(7p ? ; ~\~^ipk ff'jfe Ro ]-^»o:pCK * 0 3,^ ^ U'oc'U'ot ^ 0 

in order to satisfy realizability. The rate of change, due to the rapid-pressure correlation, of the 
eigenvalue < u a u a > is required to vanish as the limit state is approached. This insures that the 
rapid-pressure correlation model does not cause the solution to go into the unrealizable region in 
which < u a u a > is negative. This realizability limit is rephrased in terms of the determinant of 
the Reynolds stress: F = (Rjj - 3 RjjRjj + 2i?T)/6 where Rij =< U{Uj > / < u p u p > which 
can be written in terms of the invariants of the anisotropy tensor as F = 1 + 9 II + 271 1 1 where 
II = -1/2 bijbij = -1/2 < b 2 >, III = l/3bi p b p jbji = 1/3 < b 3 >. The determinant F varies 
between zero and one; F — 1 corresponds to an isotropic turbulence and F — 0 corresponds to the 
realizable limit. 
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Similar reasoning applied to the mixed tensor involving the Reynolds stress, the heat flux and the 
variance of the temperature fluctuations, produces the "joint-realizability” constraint 

D j Dt D aa ~ [^P>« ^ ][< 99 > Xj a p a < 6u a > Aj'ap] * 0 F aa ♦ 0 

which couples the rapid-pressure correlations appearing in the heat flux and the Reynolds stress 
equations. A similar determinant function Fj is defined with the normalized for which 0 < 
Fd < 1. Joint-realizability reflects the requirement that the magnitude of the correlation coefficients 
be bounded by one: the Reynolds stress and the heat-flux take on values ’’jointly” such that the 
time rate of change of D aa goes to zero as D aa goes to zero. Additional detail regarding the 
application of realizability constraints can be found in Shih and Lumley (1985). 

Note that no assumptions regarding the higher order derivatives of the eigenvalue have been made. 
The strong form of realizability, in which D 2 /Dt 2 < u a u a > > 0 is required at the realizability 
limit in order to allow the turbulence to leave the realizable state is not invoked. Such an agency is 
already present in the slow terms and it is not necessary to force the rapid terms to be responsible 
for such behavior which, as will be seen, is inconsistent with the a small parameter expansion 
of the rapid-pressure representation around the realizability limit. Instead a weak realizability 
constraint, as specified by Speciale et al. (1993), Pope (1983), which does not allow the solution 
to attain the realizable limit in finite time is invoked. This is done by requiring that the rapid- 
pressure correlation vanish more rapidly than the slow-pressure correlation model. This avoids 
any assumptions regarding the behavior of the second derivative which are required for the flow to 
leave the realizable state which is accessible in finite time in models using the strong form of the 
realizability constraint. Moreover recent work by Speziale et al. (1993) indicates that the present 
hierarchy of second-order models is inconsistent with the strong form of the realizability constraint. 
The rate of rotation of the eigenvalues coming from the second derivative is a sink term of a form 
that that cannot be balanced by the present models. In setting the portion of D/Dt < u a u a > 
due to the rapid-pressure correlation to zero while choosing a return term model that precludes 
accessibility of the realizable limit state the weak form of realizability is satisfied by the sum of the 
modeled terms on the right hand side of the transport equations. The issue raised by Speziale e< 
al. (1993) does not impact on the present model, as all the higher order derivatives vanish, and the 
realizable limit is not attainable in finite time because the relative rate of dissappearance of the 
slow-pressure terms with respect to the rapid-pressure terms. It may, however, require a rethinking 
of the modeling principles so that the intrinsic negativity of the second derivative can be properly 
balanced if the realizable limit is to be considered accessible. 

The last of the constraints to be applied to specify the constitutive relation, the geostrophic con- 
straint, is now discussed. The derivation of the constraint, as it has not been given elsewhere in 
the published literature, Ristorcelli (1991), is given. For a turbulence two-dimensional in planes 
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perpendicular to the axis of rotation Hide(1977) and Speziale (1985, 1990) have shown that the 
equations are materially-frame-indifferent. This is a direct consequence of the Taylor- Proud man 
theorem. For a modeled set of equations to be consistent with the 2DMFI principle the rapid- 
pressure correlations must satisfy the ’’geostrophic” constraints, 

^pqn^nXipqj = < UiUq > 

tpqn^nXpqi = ^ 'MqO > . 

Consider the portion of the rapid-pressure correlation associated with the rotation 

tpqk^kRo ^ X ipqj — tpqkQki^ft Ro) J U|(x)tXp(x ) > jj f q f dx /(x — X ). 

A velocity field, two-dimensional in planes perpendicular to the axis of rotation, has the representa- 
tion Up = € P qkftki>>q- Inserting the expression for the velocity field into the integral and contracting 
produces, in the integrand, the Laplacian of the streamfunction, ^, p = which reduces the 

volume integral of a two-point statistic to a local one-point statistic, 

tpqk^k R® * A ipqj — Ro * <C Ui'tp^j Ro ^tpjk^k ^ U{Up > 

upon application of Green’s theorem. The geostrophic constraint tpqn^nX ipqj — tqjn^n ^ U{Uq > 
then follows. The satisfaction of the geostrophic constraint means that, when the turbulence is 
two-dimensional in planes perpendicular to the axis of rotation, the portion of the rapid-pressure 
correlation associated with the rotation is equal to the Coriolis terms effectively removing any 
dependence on the rotation rate from the equations. Such a Coriolis-pressure force balance is 
found in the large scale atmosphere were it describes the well-known geostrophic wind. A point on 
nomenclature: a geostrophic turbulence is one that is two-dimensional because its Rossby number 
is small. This latter qualification distinguishes it from a two-dimensional turbulence (that is made 
so by some other means) for arbitrary Rossby number. Both, however, are frame-indifferent. 

The application of the five sets of constraints - normalization, continuity, realizability, joint- 
realizability and geostrophy - produce thirty-six linear algebraic equations (several of which are 
redundant) for the thirty-three unknown coefficients Ai,(i = 1,15) and D{ y (i = 1,18) appear- 
ing in the tensor polynomials. The equations are of the general form Aij(II y III)xj = 6 , where 
X{ = [Ai, ..., A 15 , ..., Dis] and where II and III are the invariants of the anisotropy tensor, 

II = -1/2 bijbij = -1/2 < b 2 >, III = 1/Sbipbpjbji = 1/3 < b 3 >. Using the definition 
F = 1 + 9 II + 27 III the general form of the equations can be rewritten as A^(//, F)xj = bi. 
Note that F and II appear linearly in the constraint equations. 

The ansatz 

Xijki = xf, kl + RX- kl 
Xijk = xf jk + FXf jk 

is used to extract more information from the constraint equations. Here X° satisfies the set of 
constraint equations A,j(//, 0 )x^ = bi obtained by application of all five sets of constraints which 
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are simultaneously valid only when F = 0 (or equivalently when III = — ( II + l/9)/3). While 
X F is obtained from the reduced set of equations Aij(II,F)xf = 6° which satisfy the three sets of 
constraints of normalization, continuity and joint-realizability and where 6° is a known function of 
the X° solution. As there are more unknowns than equations there are free parameters which will 
be used shortly to insure that the model is asymptotically consistent with an equilibrium state. 
The solution, in which all free parameters are set to zero, and which satisfy all the mathematical 
constraints, will be called the basic model and has been given in Ristorcelli (1987,1991) and Ris- 
torcelli and Lumley (1991b). These coefficients are given in Appendix 3. It is to this basic model 
that additional terms are added in order to insure that the model is asymptotically consistent with 
a known equilibrium state of a particular turbulent flow field. 

In reference to the coefficients, A t , of the basic model and defined in Appendix 3, several comments 
are appropriate. The A, are nonlinear functions of the invariants of : they are ratios of polyno- 
mials of the invariants. The realizability limit is attained along the line III = -{II + l/9)/3 on 
which F = 0. At the isotropic limit II = III = 0 and F = 1, and thus A\ + A*i = 1/10 and the well 
known exact result for isotropic turbulence is obtained. Note that the ” off-realizability” correction 
FX f is necessary to obtain this limit. The satisfaction of the isotropic limit is discovered to be 
a consequence of the constraints used to create the 2DMFI model - it is not a constraint that has 
been enforced to obtain the model but is satisfied naturally by the ansatz . This seems to vindicate 
the present procedure. 

Note that the coefficients appear to be singular at the one-dimensional limit when 1 + 37/ = 0. 
This is actually not the case as the singularities arising from the individual terms in the rapid 
pressure representation annihilate each other when summed. This has beeen verified analytically 
near the one dimensional limit, as II — ► —1/3, in principle axes of the Reynolds stress tensor using 
a perturbation expansion. 

It should be emphasized that the form of the rapid-pressure correlations was chosen so that the 
coefficients are valid for all states of the turbulence - they are not fixed to their values at the real- 
izable limit. The resulting model therefore satisfies all mathematical constraints for any arbitrary 
Reynolds stress - not just at the realizability or geostrophic limit states . The second term, FX F , 
in the expression for X, represents the ” off- realizability correction”, thus the coefficients, A{ and 
Fi, are functions dependent on the state of the turbulence as parameterized by the invariants and 
only attain their realizable limit values at F = 0. 

4. A compact representation of the 2DMFI rapid-pressure model 

The rapid-pressure correlation models are usually written in the form they appear in the Reynolds 
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stress equations with: 


IU- 2[t^ qip+£pqkftkR° l ]Xpiq 

Otj ” 2[{7^,p +6p(j/-f2/-iZo ][Xipqj -{- Xjpqi]' 

which can be rewritten in terms of the strain and rotation tensors as 

nij = 2[Sqp + Wqp][Xipqj + Xj pq i ] 

where S gp = l/2(U q ^ p + U p , g ) and W qv = [\/2{U q%p - U v>q ) + e pqk i are the usual mean strain 
rate and the total or intrinsic rotation rate. The second-order equations are then rewritten as 

D/Dt < U{Uj > +2€ t -*p < UpUj > SIRo* 1 + 2 €jk p < u p Ui > + Ylij +••• 

D/t)t < 9u{ > +2€pik$lk < 9u p > Ro ~ 1 = + Ili d" n? d”** 

where the terms omitted have already been given. Taking the contraction of the fourth-order tensor 
on the mean velocity gradients to obtain the form used in the Reynolds stress equations produces, 
with q 2 =< UpU p > y 

n-j/2q 2 = 2[A l + A 2 ]S ij 

+ [(A3 + A4 + 2As)(bi p Spj + bjpSpi ) + 2A$6ij < 65 >] 

+ [(A3 - A 4 )(b ip W pj + b JP W pi )\ 

+ [( A$ + A7 + 2Ag)(6 2 p 5 p j + 6^ p 5 P t) + 2Ag^j < 6 2 5 >] 

+ [(A 6 - M){b 2 p W p j + b 2 p Wpij\ 

+ 2[(A9 + Aio)6; P 5pg6 9 j + Ai 0 6 t j < 65 >] 

+ [(A n + A 12 + 2A 13 )(b, p S pq b 2 qj + b JP S pq b 2 qt ) + 2 A 13 (6? < bS > +b i] < b 2 S >)] 

+ [(A n - An){b ip W qp b 2 qi + b ]V W qp b 2 ,)) 

Here the angle brackets represent the trace of the indicated quantity: eg. < 65 >= 6 t; 5, ; and 
that II = — 1/2 < 6 2 >, III = 1/3 < 6 3 >. Note that the YVij has zero trace because of 
the continuity constraint, Xijjk = 0, requires A3 + A 4 + 5A 5 - //(An + A\ 2 + 4Ai 3 ) = 0 and 
Aq + A7 + 5A g + A 9 + A10 = 0 It is possible to rewrite the higher order tensor bases in terms of 
the lower order terms substantially simplifying the form of the model. The generalized Cayley- 
Hamilton theorem is used to rewrite the expression in an irreducible tensor basis. Using the matrix 
notation 

bSb = -[b 2 S + Sb 2 ]+ < bS > b + 1/2 < b 2 > S+ < b 2 S > 1 
bSb 2 + b 2 Sb = -1/3 < b 3 > S+ < b 2 S > b+ < bS > b 2 

the 2DMFI rapid-pressure correlation can be written more compactly as 

Uij /2? 2 = [B 3 + <b 2 > B>'+ < 6 3 > B?]Sii 

+ B 3 [b l p Spj + bj p S p i — 2/3 < bS > -(- B q * < bS > [b 2 j + 2///3$,j] 

+ B h [b ip W P3 + b 3P W pi ] + [B e <bS> +B 1 ' 1 < b 2 S >]b l3 
+ B 7 [b%S P3 + b 2 p S p i - 2/3 < b 2 S > 6 tj ] 

+ Bs[b 2 p W p j + bj p W p i} + B 9 [b ip W qp b 2 q} + b 3P W qp b 2 qi ] 
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where the values of the B % in terms of the Ai are given in Appendix 1. The above representation 
of the rapid- pressure correlation will be called the basic or uncalibrated model. The coefficients, 
Biy appearing in the base 2DMFI model above come from first principles: they do not result from 
any numerical optimization with experimental or numerical data. Comparisons to the FLT model, 
Fu et al (1987), shows that the two models have the same tensor structure containing the same 
generators. The only difference is in the terms proportional to [6*j] and \b 2 - + 2///3^j], usually 
identified with the slow-pressure’s contribution to the pressure-strain correlation and calibrated 
accordingly, now reflect a contribution from the rapid-pressure. Contraction of the irreducible form 
of Xijki on the mean velocity gradients has produced the bases [b{j] and [b 2 j + 2JJ/3£,j] in which 
the coefficients are functions of the invariants < bS > and < b 2 S > in addition to the dependence 
on < b 2 > and < 6 3 > appearing in the A{. The two tensor bases, usually associated with the 
slow-pressure correlation, arise as a consequence of starting with the two-point volume integral 
and represent a contribution to the pressure-strain correlation whose structure is identical to the 
slow-pressure models but whose genesis is in the rapid-pressure. 

Speziale et al. (1991) have also written a general form for the pressure-strain covariance. It is linear 
in the mean velocity gradients and nonlinear in the anisotropy tensor. Their expression contains 
the same generators as the present model except for the cubic term [bWb 2 — b 2 Wb]. Speziale 
et al (1991) use results from rational mechanics (cf. Smith (1971)) to expand in a functional 
basis. The present strategy uses a polynomial basis for which the generator [bWb 2 — b 2 Wb] is 
not redundant, Spencer (1971). The tensor polynomial given above is irreducible and the basis is 
optimal. Speziale et al (1991) have shown, for planar flows, that the generators nonlinear in the 
anisotropy tensor can be expressed in terms of generators linear in the anisotropy tensor. This fact 
led to the very simple form of the SSG model, Speziale et al (1991). This fact, from a rigorous 
though not necessarily practical point of view, limits the model to planar flows. The results of 
Speziale et al (1991) can be used to recast the present model into its linear planar form. For 
planar flows Speziale et al (1991) has shown that 

[bfpSpj 4- bj p S p i - 2/3 < b 2 S > <5 t ,] = -b 33 {b ip S PJ 4 b jp S p i - 2/3 < bS > - 2j3(III/b 33 )S,j 

[b] p W PJ 4 bjpWpi] = -b 33 (b, p W pj + b JP W pi ] 

from which it follows that 

[b tp W qp b 2 qj 4 b }P W qp b 2 qt ) = (II 4 b 33 b 33 )[b tp W P} + bj p W p i] 

and the planar form of the 2DMFI rapid-pressure model can be written as 

nr, /2 g 2 = (C 3 - 2IIC 3 + 3/7/Cr - 2/3(III/b 33 )Cj)(S t] \ 

+ (C4 — b33C7)[bi p Spj + b 3 pS v% — 2/3 < bS > 6 t y] + C™ < bS > [tifj + 2///3<5 t y] 

+ (C5 — ^33^8 + ( II + ^33 ^33 ) C9 ) [6|p fF pj + bjpWpi] + (^6 < bS > +CT < b 2 S >)[f>ij] 
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which has the same linear tensor bases as the quasi-linear SSG model. The form of the model 
in planar flows is consistent with to the SSG model which is the topologically generic form of a 
general class of models for the fourth order rapid-pressure tensor X{jk\ . The coefficients in the 
model, however, are nonlinear functions of the invariants while in the SSG model, given below for 
comparison, 

TVij = — (2Ci£ + ClV^bii) + C 2 t\b\ + 211/Mij] + (C 3 - (-2 :JJ)»/ac3)*[Sy] 

+ C 4 fc[6, p 5pj + bjpS p i — 2/3 < bS > ft tj ] + C$k\bipWpj + bjpWpi] 

the Ci are numerical constants and one coefficient in the [Sij] term, in the rapid-pressure portion of 
the model is nonlinear. The C, in the SSG model are determined by a numerical optimization so that 
the model reproduces as closely as possible: 1) the stationary state of the homogeneous shear and 
2) maximizes the kinetic energy growth rate of the rotating homogeneous shear as close as possible 
to the Sl/S = 0.25 predicted by rapid distortion theory without introducing a Richardson number 
similarity, Speziale and Mhuiris (1989), while insuring that the points of exchange of stability are 
outside of those predicted by the linear theory of Bertoglio (1982). 

In the return term of the SSG model the correction to the linear term arrived at more or 
less intuitively by Speziale et al. (1991) is, to lowest order, vindicated by the present results. 
Speziale et al. (1991) have altered the term linear in 6;j, usually associated with the slow-pressure 
correlation, to include a term involving the mean flow, a term proportional to the production 
V = -2 < UiUj > Ui,j\ the usual C\sbi 3 — ► (C\S + C{V)bij. The present analysis indicates that the 
portion of the rapid-pressure contribution to the pressure-strain correlation, linear in b{j y has the 
form [Bq < bS > +B^ f < b 2 S >]6,y. Note that < bS > can be written in terms of the production 
as < bS >= — V !q 2 . The present analysis suggests the possibility of adjusting the nonlinear return 
term for mean velocity gradient effects in a similar way. In the present model the portion of the 
rapid-pressure quadratic in the anisotropy tensor is B f ± <bS > [tifj -I- 2/7/36, y] and therefore also 
scales with the production. 

Speziale et al. (1992) have reflected on the ambiguity of the distinction between the rapid- and slow- 
pressure contributions to the pressure-strain correlation. They have modeled the whole pressure- 
strain correlation but whether to interpret their adjustments (as described in the previous para- 
graph) as incorporating the effects of the mean strain on the return terms or as contributions of the 
rapid-pressure to the total pressure-strain is not clear. Here, however, the distinction is clear. The 
present analysis starts from the tensor polynomial representation for the rapid-pressure integral 
and produces terms whose tensor structure is identical to that which are traditionally called the 
slow-pressure. The other non-linear pressure-strain models, FLT or the SL model, Shih and Lumley 
(1985), do not have terms in or b 2 - that can be similarly identified. 

The derivation of a rapid-pressure correlation satisfying all the limit states and valid away from the 
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limit state, though mathematically rigorous, will not necessarily produce a model that performs 
better than other models in flows for which latter models are calibrated. However, for the general 
class of flows for which second-order models are suitable we are of the opinion that the satisfaction 
of all mathematical constraints constitutes a necessary (though not sufficient) requirement to assure 
the predictive capabilities in flows different from those for which the models are calibrated. The 
existence of several free parameters can then be used to calibrate the model to specific classes of 
flows of computational interest. This tack is taken in the next section where the fixed points of the 
modeled equation are matched to the fixed points of the homogeneous shear. In this way a model 
for the class of flows in which the mean shear is the dominant production mechanism is created. 

5, Calibrating the rapid-pressure correlation representation 

The coefficients, appearing in the 2DMFI model come from first principles: they do not result 
from any calibration with experimental or numerical data. They represent the minimum number 
of determined coefficients necessary to satisfy the mathematical constraints on the rapid-pressure 
correlation for a three-dimensional Reynolds stress. It, however, is not a unique representation: 
there are an infinity of solutions corresponding to different values of the free parameters which, in 
the basic model representation shown, have been set to zero. Computations have shown that the 
predictive capabilities of the basic model are inadequate. The mathematics built into the model do 
not capture the experimentally known stationary points of homogeneous shear. To compensate for 
some of the approximations made in the mathematical development the model is modified so that 
the modeled evolution equations have the same asymptotic behavior as that observed: to the basic 
model additional terms are added to insure that the model is consistent with an equilibrium state 
of a particular benchmark turbulent field. This is done without sacrificing any of the mathematical 
principles built into the model. It, however, can not be done arbitrarily. 

The strategy is to require asymptotic consistency with an equilibrium state. The modeled evolu- 
tion equations are required to have the same fixed points as those observed in experiment. The 
equilibrium constraint is used to obtain additional constraints equations to specify the free param- 
eters. This is very similar to the strategy employed so far in that limit states are used to set model 
coefficients. Here, of course, the equilibrium state is much closer to those expected to be seen in 
flows of engineering interest. The free parameters will now be called calibration coefficients, Af, 
and will be collectively denoted by X?? k i appearing in the decomposition 

Xijki = Xf- kl 4 FX[- kl 4* FX-f kl . 

X°? k i represents the additional constant terms necessary to capture the stationary state. This form 
is equivalent to assuming that the coefficients in the tensor polynomial have the form A, = A® + 
FAf 4 FA^. The A\ satisfy normalization and continuity constraints and the algebraic equations 
describing the stationary state of a homogeneous shear. It is at this point that numerical or 
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experimental data for the asymptotic values of the anisotropy tensor, the production to dissipation, 
(Pie) oo, and the ratio of time scales, (5 , fc/e) 00 , are necessary. The details of this strategy are now 
given. 

The A£ satisfy the six equations given by the homogeneous form of normalization and continuity 
constraints, X^^ — 0, — 0, which are valid for all states of the turbulence. Substituting in 

the calibration coefficients allows six .4 * , to be expressed in terms of seven free parameters, 

A\ = -(50Ag + 12Ag + 4i4j 0 )///15 + ( A\ x + A c n - 6A 13 )/77/5 
A% = (20 Ag + 3Ag + A C 10 )II/15 - (3 A c u + 3 A\ 2 + 2A 13 )/77/10 
Ag = -11/3A| + (A c n +3AU + 8A\ 3 )II/3 
Ag = -4/3 Ag + 2 (A c n + 2Ag 3 )///3 
A§ = -l/3(llAg + 3Ag + Ag 0 ) 

A$= — 2/3(2Ag + i4j 0 ). 

The seven parameters will be determined by asymptotic consistency with a particular equilibrium 
state. At this point the rapid-pressure correlation model is fully general and it is possible to write 
it in its final form, without specifying the calibration coefficients, as 

mi/29 2 = [C 3 - 2//C3 + 3IIIC" r )Sij 

+ C 4 [b xp s pj + b 3P S pi - 2/3 <bS> Sij] + C'" < bS > [&?■ + 277/3<5 t j] 

+ C 5 [b xp w pj + bj p W pt ] + [Ce <bs> +Cf < b 2 S >]b tJ 
+ C 7 [b%S pj + bj p S p i - 2/3 < b 2 S > S i3 ] 

+ C s [bf p W pj + b%W pi ] + C 9 [b ip W qp b 2 q] + b Jp W qp b%). 

The calibration coefficients, A g, have been added to the base model coefficients, 77,, to obtain the 
final model coefficients C,: 

C 3 = B 3 - 2F(10Ag + 3 Ag + Ag 0 )///5 - F(A C U + A c u + 14Ag 3 ))77//5 
C'i= B'i + F(A% + A\ 0 ) 

B'l'- \!3F(A\ X + Ag 2 + 2Ag 3 ) 

C a = Ba + F(-3A\ + II(A‘ n + A\ 2 + 4Ag 3 )) 

C"' = B'Z + F(A c n + Ag 2 + 4 Ay 
C 5 = Bs + F(-7/3Ag + ( -A c n + 3 Ag 2 + 4Ag 3 )///3) 

C 6 = fl 6 + F(2Ag + 4Ag 0 ) 

C 7 = B 7 - 3F(Ag + Ag + Afo) 

C» = B s - l/3F(7Ag + 3Ag - Ag 0 ) 

C 9 = B 9 + F(A c n - A\ 2 ) 

Note that the traces can be written in terms of the production < bS >= —V/q 2 and < b 2 S > = 
1/2 b 33 P/q 2 ', in planar flow < bS >= — 2 b\ 2 S and < b 2 S >= —b\ 2 b 33 V /q 2 . 

The calibration coefficient Ag involves adjustments to the generators [bS + Sb - 2/3 < bS > 1] 
and [bW — Wb] which are linear in the anisotropy tensor. Numerical experiments have indicated 
that it is important in establishing the levels of the normal components of the Reynolds stresses. 
Combinations of Ag, A% and Ag 0 affect the generators quadratic in the anisotropy tensor, [b 2 S + 
b 2 S — 2/3 < b 2 S > 1] and [b 2 W - Wb 2 ] and to a very small degree the S term. Combinations of 
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only Ag and A xo can be used to control the contribution to Y\ r proportional to b and S. Experience 
with the FLT model, Fu et ai (1987) indicates that the cubic term [bWb 2 - b 2 Wb] involving 
the rotation tensor is important in controlling the relative level of the normal components of the 
Reynolds stress in situations with rotation. This suggests that the calibration coefficients A\ X ,A\ 2 
and A \ 3 are important. The combination of calibration coefficients A\ x — A\ 2 controls tensor 
products like [bWb 2 — b 2 Wb] and combinations of terms A\ X ,A\ 2 and A X3 control the contributions 
of the b 2 and S. 

The existence of the free parameters allows the model to be calibrated to a specific class of flows 
of computational interest. For example, in a buoyant flow one might evaluate the fixed points 
using the experiments of the homogeneous shear in a constant mean temperature gradient thus 
producing a model suitable for a class of stratified flows of geophysical interest. This calibration 
would, of course, involve the model for the rapid-pressure correlation appearing in the heat flux 
equation. For many engineering problems the primary production mechanism is the mean shear 
and capturing the fixed points of the homogeneous shear in the modeled equations will make the 
model suitable for a wide class of flows. At this point one could also consider using the exact 
results of rapid distortion theory to obtain values for the calibration coefficients. This, however, is 
inconsistent with the equilibrium hypothesis underlying the local approximation to the constitutive 
relation invoked for the rapid-pressure correlation integral. Moreover the rapid distortion problem 
is linear admitting a superposability that is not possible in the context of nonlinear second-order 
closures. The homogeneous shear structural equilibrium will be used to fix the representation for 
the rapid-pressure correlation. 

The fixed points of the homogeneous shear are now built into the model by specifying the calibration 
coefficients. The modeled evolution equations for the homogeneous turbulence in a mean velocity 
gradient, are 

D j Dt <c 'Ui'iij — 2 tikp *C UpUj 2 tjkp ^ ^p^t ^ ^ u^Up > JJ j ,p <c ujitp U%^p 

+ nJi - C x eb X] + C 2 e[bl } + 2II6ij/3] - 2/3e6ij 

D/Dt e = -(C £ i < UiUj > Ui,j + C e 2£) e/k 

assuming local isotropy for the dissipation. The terms -C\£b{j + C 2 £[bjj + 2/7/36,y] represent the 
return to isotropy pressure correlation. Using < u t Uj >= q 2 (b{j + 1/3 the equations for the 
anisotropy are 

D / Dt b tJ — Qtjkpbpi [bipUj,p ~\~b 3 pUi,p 2/36, y < bS >] 

—2/3 Sij + 26, j < 65 > + n r ij /q 2 ~ ( C , - 2)6, ye/? 2 + C 2 [6?- + 2/7/36, y] £ / 9 2 . 

The mean strain and rotation tensors are 5,y = l/25*(6,i6y 2 + 6, 2 6yi) and W tJ = l/2W*(6n6j2 — 
6, 2 6yj), where, 5 = f/i, 2 = S* = W*. Setting the D/Dt bij = 0 produces three algebraic equations. 
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Inserting the experimentally determined asymptotic values of bij and Sk/e produces three addi- 
tional constraints leaving four of the seven A * as free parameters. The data from the experiments 
of Tavoularis and Corrsin (1981), Champagne, Harris and Corrsin (1970), Tavoularis and Karnik 
(1989), and the DNS of Rogers et ai (1986) are summarized in the Table 1. 



TC 

CHC 

TK.A 

TK.C 

TK.D 

TK.G 

TK.J 

TK.K 

DNS 

bn 

0.197 

0.137 

0.217 

0.257 

0.197 

0.157 

0.157 

0.147 

0.215 

bn 1 

-0.14 

-0.165 

-0.165 

-0.165 

-0.17 

-0.148 

-0.154 

-0.149 

-0.158 

b 22 

-0.143 

-0.083 

-0.133 

-0.143 

-0.133 

-0.113 

-0.103 

-0.093 

-0.153 

^33 

-0.053 

-0.053 

-0.083 

-0.113 

-0.063 

-0.043 

-0.053 

-0.053 

-0.062 

(P/e) 00 

1.75 

1.0 

1.38 

1.37 

1.64 

1.33 

1.45 

1.37 

1.80 

{SK/e) 00 

6.25 

3.03 

4.2 

4.15 

4.82 

4.5 

4.71 

4.60 

5.7 


Table 1. Data for the homogeneous shear flow. 

There is considerable scatter in the data due to the technique used to generate the turbulence, 
individual wind tunnels in which the different experiments were done, and experimental error. Not 
all the data, as has been qualified in the references from which the data is drawn, represent the 
asymptotic state. The following fixed point values are taken to be representative: b^ = 0.203, 
b? 2 = -0.156, bf 2 = -0.143, b? 3 = -0.06, (Sk/e) oo = 5.54, (P/e) <x> = 1.73. They are obtained by a 
simple average of the data of TC, TK.D and the DNS. These three cases are chosen because they 
have the highest values of of the nondimensional time (SK/e) <*> corresponding to the flows that are 
furthest in their development to the asymptotic state. The values of the invariants corresponding to 
these values of the anisotropy tensor are: //oo — — 0.057, ///oo — 0.0032, Foo = 0.573. Substituting 
these asymptotic values into the fixed point equations, D/Dt b{j — 0, reduces the number of free 
parameters from seven to four. 

A% = -0.29 - 0.06(^ 0 - A c s ) 

A c u = -3.6 + 5A C 10 - 2 A c 13 - 12.7i4g - 3.8Ag 
A\ 2 = -24.5 - 44.2i4f 0 - 2 A\ z + 29Ag - 8 A% 

This set of constraint equations is dependent on the model for the return to isotropy pressure. For 
simplicity the the nonlinear return coefficient has been set to zero, C 2 = 0, and the well accepted 
value, C%° = 3.4, has been chosen. There are still four free parameters. The following values of the 
four free parameters are chosen: Ag = 0.8, Ag = — 1.0, AJq = 0.01, A\ 3 = 0. The undetermined free 
parameters have been set by matching the values of the anisotropy for the log layer. The procedure 
is outlined in more detail in Appendix 1. 

6. Computations and comparisons in homogeneous turbulence 

The 2DMFI model falls into the same class of representations as the FLT and SL models: they all 
use nonlinear terms and invoke some form of realizability constraint to evaluate the coefficients. For 
this reason the 2DMFI model will be compared primarily to the nonlinear SL and FLT models. For 


17 



completeness and because it appears to be a very successful model for the planar flows, computations 
with the quasi-linear SSG model are also shown. The SSG model is linear in the anisotropy tensor 
though nonlinear in that the scalar coefficients are functions of the invariants of the the anisotropy 
tensor. It should, however, be kept in mind that the SSG model satisfies realizability for the kinetic 
energy and not for the individual Reynolds stresses and is therefore in another class of models. This 
issue is more fully explored in a subsequent section. Results are not compared to the LRR model 
as, in concordance with the observations of Speziale et al (1991), the SSG model is viewed as an 
updated optimized LRR model. 

In all the calculations with the 2DMFI model a simple linear Rotta type model for the slow pressure 
correlation will be used. This corresponds to C 2 = 0 in the canonical form given above and is 
consistent with the present calibration to the homogeneous shear. For the linear return coefficient 
a simple expression, C\ — 2 — 31 IIF 1 / 2 , is used. This satisfies the isotropic limit, C\ = 2.0, and 
is consistent with the assumed value for the asymptotic homogeneous shear, = 3.4. The form 
chosen is consistent with a weak form of realizability and the recent results of Speziale et al, (1993) 
regarding the rate of disappearance of the rapid-pressure correlation relative to the return pressure 
correlation as the realizability limit is approached. 

The values used for the constants in the dissipation equation are: C € \ = 1.44, C e 2 = 1.83. Note 
that this corresponds to a single universal fixed point (V/e) <*> = 1.88 independent of rotation. 
The single fixed point is a well-known deficiency common to all the present forms of the modeled 
dissipation equations. 

Case 1: Homogeneous shear 

The calibrated model is now used to compute the time evolution of the homogeneous shear flow. The 
mean strain and rotation tensors are Sij = l/2S*(Mj 2 + Mji) and Wij = l/2W*(6 tl 6 j2 - S i2 Sji), 
where, S = U\, 2 = S* = W* . In Figure 1 the time evolution of the turbulence energy is compared 
to the LES of Bardina et al (1983), and the three models FLT, SL and SSG. A similar monotonic 
behavior is found for other statistics, b\ 2} &n,// in the flow and, as they do not constitute new or 
different information, are not shown. In general, starting from physically realistic initial conditions, 
the flow attains its asymptotic state rapidly and monotonically. The asymptotic states which the 
different models attain are given in the accompanying table. The column labeled experimental data 
is an average of the three cases TC, DNS and TK.D. 
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Equilibrium 

Values 

2DMFI 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

Experimental 

Data 

LO° 

°n 

0.209 

0.202 

0.208 

0.219 

0.203 

h°° 

°12 

-0.155 

-0.080 

-0.146 

-0.164 

-0.156 

h°° 

°22 

-0.148 

-0.195 

-0.144 

-0.146 

-0.143 


-0.061 

0.007 

-0.064 

-0.073 

-0.06 

(P/e) 00 

1.88 

3.42 

1.99 

1.88 

1.73 

(SK/e) 00 

6.08 

21.35 

6.84 

5.76 

5.54 


Table 2. Comparison of the model predictions for the equilibrium values in homogeneous shear flow 
(V/e = 1.88) with the experimental data given in Tavoularis and Karnik (1989). 

Note that the different ( V/e )oo attained are functions of the different C c \ and C z 2 used in the models. 
The present form of the dissipation equation insures that the quantity (V/e )# 0 = 1-88 regardless 
of initial conditions for all Sk/e. This is a shortcoming of the modeled dissipation equation and 
shows up in a larger bn than the flow for which it was calibrated in which (V/e)^ = 1.73 

Case 2: The equilibrium wall layer 

Another simple but important test cases is whether the model can capture the stationary state of 
the log-layer in the channel flow. The homogeneous shear and the log-layer are similar in that they 
achieve, to a suitable approximation, an equilibrium state. Abid and Speziale (1993) have discussed 
the relevance of this test case and noted the inability of most rapid-pressure closures to perform 
successfully in the log-layer. Our results are in agreement with their contention that a model which 
is asymptotically consistent with the stationary states of the homogeneous shear will also do well 
in the log-layer. The models are compared to the DNS of the channel flow of Kim (1993) which 
is an update of the simulations reported in Kim et al. (1987). The data presented represent an 
average of the values of the anisotropy in the region 70 < < 100 outside the viscous sublayer. 


Equilibrium 

Values 

2DMFI 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

DNS 

Data 

Experimental 

Data 


0.180 

0.079 

0.141 

0.201 

0.180 

0.22 


-0.141 

-0.116 

-0.162 

-0.160 

-0.134 

-0.16 

LOO 

°22 

-0.142 

-0.082 

-0.099 

-0.127 

-0.140 

-0.143 

^33 

-0.039 

0.003 

-0.042 

-0.074 

-0.040 

-0.06 

(P/e) 00 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

(SK/e) 00 

3.55 

4.30 

3.09 

3.12 

3.73 

3.1 


Table 3. Comparison of the model predictions for the equilibrium values in the log-layer of turbulent 
channel flow (V/e = 1) with the DNS data of Kim (1993) and the data of Laufer (1951) given in 
Abid and Speziale (1993). 

Case 3: Homogeneous shear with rotation 

The present test case, the homogeneous shear with rotation, and the next test case, the homogeneous 
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shear with streamline curvature, are important test cases. In both these cases additional forces, 
which stabilize or destabilize the flow, are present. These effects appear in the evolution equations as 
additional production mechanisms for the Reynolds stresses. In the case of rotation the production 
terms in the evolution equation for the turbulence kinetic energy do not directly depend on the 
rotation: the turbulence energy production depends on the rotation only through the off-diagonal 
components of the Reynolds stress while in the case of streamline curvature the production terms 
in the evolution equation of the turbulence kinetic energy do directly dependent on the curvature. 
These two cases are important test cases not only because 1) the models have not been calibrated for 
them but also because 2) the model will have to predict both the stabilization and the destabilization 
of the turbulence and 3) the critical values of the governing parameters which demarcate the regions 
of flow stabilization from flow destabilization. 

For flows in the rotating coordinate system the Coriolis terms must be carried and W{j appearing 
in must be replaced by the total rotation tensor Wij + ejik Cl^. Thus in the rapid-pressure model 

W* = 5(1-211/5). 

Figure 2b-2e show how the models perform in rotating shear for rotation to shear ratios Cl/S = 
0,1/4, 1/2 compared to the LES data of Bardina et al . (1983) in Figure 2a. In general, all the 
models are able to capture both flow stabilization for some A > Cl/ S > B and flow destabilization 
for some A < Cl/S < B. The points A and B represent the points of neutral stability on a 
bifurcation diagram in the phase plane (e/Sk) oo and (Cl/S)^. All of the models have a bifurcation 
diagram of the same general form, Speziale and Mhuiris (1989), indicating a stabilization of the 
flow outside of some region of approximate size 0 < Cl/S < 0.5, predicted by the linear rapid 
distortion theory of Bertoglio (1982). The most important facts concerning the different models 
for the homogeneous rotating shear can be summarized by indicating the unstable regions in which 
the models predict a non- trivial equilibrium ( e/Sk ) co . 


SSG: 

-0.09 < 

Ll/S < 0.53 

RDT : 

0.00 < 

a/s < o.5o 

2DMFI : 

-0.063 < 

a/S < 0.502 

SL : 

-0.14 < 

a/s < o.4o 

FLT : 

-0.11 < 

a/S < 0.39 


Near the point of linear neutral stability Cl/S = 1/2 both the SL and the FLT models predict a 
premature restabilization at values of Cl/S 20% and 22% lower than that predicted by the linear 
theory. The 2DMFI model is within 4% of the linear prediction. 

None of the models tested, linear or nonlinear, captures the point of maximum kinetic energy 
growth at ( Cl/S) ma x = 0.25. To do so would mean that the equations would exhibit a Richardson 
number similarity which, as Speziale and Mhuiris (1989) have shown, is not admitted by the Navier 
Stokes equations. The two models that come closest to (Cl/S) max — 0.25 are SSG at Cl/S — 0.22, 
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which was calibrated using this fact, and 2DMFI at fi/5 = 0.2 which was not calibrated using 
any rotating flows. The current modeled dissipation equations predict a (V/e)co = ~2b\2Skje = 
(C e 2 — l)/(C e i — 1) = const., where the constant is model dependent but independent of rotation 
rate. The constant attains the value for the equilibrium homogeneous shear for arbitrary rotation 
rate, a fact which is not consistent with observation. In a flow that is stabilized by rotation, say 
11/5 = 1, production must be less than dissipation for the equilibrium state to be reached. The 
dissipation equation cannot be used for calibration in rotating flows without compromising the 
model when a dissipation equation capable of predicting the stationary values of ( Vje ) QO = f(fl/S) 
becomes available. 

Case 4: Homogeneous shear with curvature 

For the homogeneous shear with streamline curvature the mean strain and rotation tensors are 5,j = 
l/25*(Mj 2 +Wji) and W {j = l/2W*(6n6 j2 - S i2 S jx ) where 5* = S(l-atb) and W* = 5(1 +stb) 
where stb = ( U c /R c )/S is the stability parameter. The geometry for the curved homogeneous shear 
follows that of Holloway and Tavoularis (1992): R c is the radius of curvature of the flow and U c 
is the axial velocity at the centerline and the crosstream gradient of the axial velocity is the shear 
Uiy 2 = S' The stability parameter has been renamed stb so as not to confuse it with 5 which is 
traditionally used for the mean shear. The kinetic energy growth rate is suppressed, relative to the 
homogeneous shear, for stb > 0 and increased for stb < 0, while for stb > 0.05 the experimental 
data indicates a relaminarization. 

Figure 3 compares the model results to the experimental data for 612 versus stb. The plot has been 
generated by computing the flow from the beginning of the straight section of the wind tunnel to a 
value of St = 10 which corresponds to the end of the curved section. The initial conditions on the 
second order moments are given by the experimental data. The initial condition on the dissipation 
rate is determined by matching to the kinetic energy growth rate at the beginning of the straight 
section. 

The different models all capture the trend in the stabilization/destabilization with respect to the 
stability parameter. The primary difference in the predictions of the different models seems related 
to their ability to capture the homogeneous shear at stb = 0. The results of the different models 
would be in more agreement for negative and small positive values of stb if they predicted the same 
results for the homogeneous shear. 

The bifurcation diagram for the second order models, in Figure 4, was generated by letting the 
solution procedure go to its asymptotic state. There is a critical value of stb c at which the stabi- 
lizing effects of curvature begins to causes a negative kinetic energy growth rate, which ultimately 
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relaminarizes the flow. The critical values, stb c , predicted by the different models, are 


HkT : stb c = 0.05 

2 DM FI : stb c = 0.067 

FLT : = 0.075 

SSG : = 0.10 

5L : = 0.105 


where HkT is from the experimental data of Holloway and Tavoularis (1991). There is a consistent 
trend for the SSG, 2DMFI and FLT models, when compared to the critical values for the rotating 
shear: the higher (f 1/5% for stabilization of the flow correspond to higher stb c . The SL model has 
small but nonzero e/Sk over the range 0.25 < stb < $tb c . 

Case 5: Two and three-dimensional strains 

The rapid pressure model is used to compute three strain flows: the plane strain and the axisym- 
metric contraction and expansion. These flows are another test case as the rapid pressure model 
has not used these flows to set the calibration coefficients. The results are compared to the direct 
numerical simulations of Lee and Reynolds (1985). Because the simulations are conducted at low 
Reynolds number the anisotropies are expected to be somewhat higher than those of a fully de- 
veloped turbuelence. However the use of the physical experiments conducted at higher Reynolds 
number is also somewhat tenuous as the initial conditions on e/Sfc, as has been pointed out by 
Speziale et al. (1991), are not known with certainty. The same test cases as those given in Speziale 
et al. (1991) are used. Our results are also compared to the SSG model as it appears to be the 
current model that gives the best results. The evolution of the kinetic energy for these flows is not 
presented; the results for the models are in very good agreement with the data and each other and 
do not definitively distinguish between the various models. 

Figure 5 shows the evolution of the anisotropy for the plane strain, = S*(6n8ji — ^ 2 ^ 2 ) 
starting from isotropic initial conditions. Figure 6 and 7 shows the evolution of the anisotropy for 
the axisymmetric contraction and expansion. Here for the contraction 5^* = S*(6n6ji - l/2^2^2 “ 
1/2 6 i3 6 j3 ). For the expansion S is replacesd with — S . Results for ail the plane and axisymmetric 
strain flows capture the trends nicely. 

8. General discussion 

The constraints of geostrophy, realizability, joint-realizability, normalization and continuity have 
been used to create the 2DMFI model. The use of a realizability type constraint to obtain values 
for unknown coefficients in the models, in the case of the strong form of realizability, has been 
criticized on the grounds that one should not use an extreme state to set the coefficients in an 
equilibrium model. This certainly is the case in any rigorous interpretation of the statement - 
requiring all the scales of the motion to satisfy the indicated limit. However, from the point of 
view of a useful engineering approximation, it should be kept in mind that < UiUj > represents 
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an integral over all scales of the motion: from the production scales at kI 1 and larger to the 
dissipation scales kt} ~ 1. In a turbulence with a k“ 5 / 3 inertial subrange in which there is enough 
of a separation of scales for a second-order simulation to be useful, say at least Rei ~ 10 4 , the ratio 

Q / A 

between the dissipative and the energy containing length scales is rj/l ~ Re t ~ 1000 and the 
flow scales range over 0 < kI < 1000. However, approximately 85% of the energy of the motion is 
contained in the first decade kI < 10: the major contribution to < U{Uj > is from the scales of the 
motion greater than one tenth of the production scales, for which k^I ~ 1 - If only the largest 1% 
of the flow scales, i.e. from 0 < ni < 10, begins to lose an eigenvalue of the Reynolds stress tensor, 
through some dynamical or kinematical agency, < U{Uj > begins to approach the realizable limit. 

Part of the reluctance to accept the second-order modeling technique is that turbulence simulations 
of complex inhomogeneous flows with multi-dimensional mean flows which may have body forces, 
streamline curvature or rotational effects, are extremely difficult to compute. ^From a strictly 
practical point of view - computability - incorporating the realizability constraints into the models 
for unknown correlations has some very tangible and beneficial effects. During the convergence to 
a solution, from more or less arbitrary initial conditions, the iteration will be plagued with negative 
normal stresses and correlation coefficients larger than one. When this occurrs the solution is clipped 
and the solution procedure restarted from the new clipped initial conditions. The frequency of this 
clipping and resetting procedure is substantially reduced (in simple flows even eliminated) when 
using realizable models. 

How one gets to a realizable turbulence is not an issue in problems with steady states, as long as 
the final state is realizable. This, however, is not the case for problems that are unsteady. For time- 
varying flows, varying on the integral time-scale, the quasi-steady problem for which the second- 
order methods are still suitable - this issue becomes serious. Our experience with the Reynolds 
averaging procedure, in buoyantly driven elliptic flows with rotation, indicates that it acts as a low- 
pass filter: the rapidly fluctuating instantaneous dynamics are subsumed by the averaging procedure 
leaving the slow-time large-scale parts of the flow, evolving on time scales commmensurate with the 
integral time scale, to be captured by the computation. If the simulation is to reflect the physics 
of the time evolution of the flow, it must stay realizable. Clearly excessive realizability violations, 
requiring a clipping and resetting of the solution, which produces a solution that never evolves far 
from the transient associated with the most recent clipped initial condition, are not acceptable. 
In such flows, satisfying the realizability constraint has very important consequences regarding the 
validity of the time evolution of the flow. 

The coefficients of the tensor polynomials used to represent an unknown correlation in a constitutive 
relation are, according to the theory, non-constant functions of the invariants of the independent 
tensors and thus depend on the state of the turbulence. These coefficients, in ’’realizable” turbulence 
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models using some form of the strong realizability constraint, are obtained by requiring that the 
rate of change of a positive semi-definite quantity be zero when some limit state is achieved. The 
coefficients so obtained are constants and are strictly valid only at the limit state. To the limit state 
value of the coefficients, which are regarded as a zeroeth order approximation to the coefficient in 
a general flow, ad hoc corrections are then added. These corrections, which vanish as the limit 
state is approached, require some sort of numerical optimization with well known simple flows. 
This summarizes the methodology used in other ’’realizable” turbulence models. In the linear SSG 
model the coefficients are, for the most part, constants that are also set by matching to a limit state. 
In the case of the SSG model the structural equilibrium limit is used to set the constants in the 
model. This is a much closer approximation to the turbulence expected to be seen in engineering 
problems. The constants in the SSG model may be viewed as the values the nonconstant coefficients 
have near the equilibrium state. 

In the present method the realizability principles are used to obtain non-constant values of the 
coefficients in the constitutive relations that are valid for all states of the mechanical turbulence. 
Recall that the basic form of the rapid-pressure model is comprised of two parts Xfj k( and FXf jkl , 
the ’’off- realizability” correction. X® jkl satisfies simultaneously the five constraints - geostrophy, 
realizability, joint-realizability, normalization and continuity, while FX^ kl , also obtained analyt- 
ically, satisfies the substantially less extreme joint-realizability as well as the homogeneous form 
of the normalization and continuity constraints. Thus, although the model is consistent with the 
realizable limit state, the values of the model coefficients are not the values the coefficients have at 
the limit state. An extreme state of the flow has only been used to set the coefficients in the Xf jkl 
part of the model - the additional FX^ kl terms are, for the mechanical turbulence, fully general. 

It is to this basic model, valid for all states of turbulence not just at the realizability limit, that 
one adds the FXf? ki term that is necessary to obtain the equilibrium values in the flows evolution 
from arbitrary initial conditions. The requirement of asymptotic consistency with an equilibrium 
state, first used by Speziale et al (1991), is the single most important and consistent physical 
requirement one can impose. Second-order closure methodology is built around the assumption 
that there is, to a suitable approximation, for the class of flows to which second-order methods 
are appropriate, a equilibrium state and in the absence of disturbing forces the flow relaxes to 
that state on a time scale similar to the eddy turnover time. It is precisely this phenomenological 
behavior that is built into the model by requiring the fixed points of the modeled equations to be 
consistent with those obtained from experiment. The assumption that allows the parameterization 
of the two-point correlation as a local function of the anisotropy tensor, XijM = Xijkl(bij), Is such 
an equilibrium assumption. 

The penalty paid for these additional features associated with the satisfaction of the mathematical 


24 



constraints is a more complex model. It should, however, be pointed out that the present 2DMFI 
model has the same tensor basis as the FLT model and is therefore no more complex except for 
the expressions for the nonconstant coefficients. Moreover, the penalty is slight in the light of 
the reduction of the computational difficulties found during the calculation of quasi-steady time 
evolving flows with this representation for the rapid-pressure. The model, along with several other 
models, has been used to compute inhomogeneous buoyancy driven rotating flows that occur in 
the Czochralski crystal growth melt in which the Reynolds stress are three-dimensional, Ristorcelli 
and Lumley (1991a, 1993), Ristorcelli (1991). In computing these time varying flows it was found 
that the present 2DMFI model produced far fewer realizability violations during the course of the 
flows evolution. For a quasi-steady flow this is a crucial point: every time realizability is violated 
the solution is reset and the solution never evolves past the transients associated with reseting the 
initial conditions. Such a computation cannot be expected to reflect an ensemble average of the 
original system. 

9. Suggestions for future work 

In the effort to produce a representation for the rapid-pressure correlation valid for the class of 
flows to which second-order modeling is suitable, some shortcomings in the data on homogeneous 
’’building block” flows have become apparent. Though the homogeneous shear seems reasonably 
well-documented it is not clear that the asymptotic states have been reached in some of the exper- 
iments. Moreover the discrepancy between the high values of bfi obtained in the DNS versus those 
seen in the laboratory data has not been explained. Additional work expanding on the notion of 
two classes of flows, according to Tavoularis and Karnik (1989), and therefore the possibility of two 
equilibrium states, merits investigation. 

For the homogeneous shear with rotation, a very basic flow, there doesn’t seem to be substantial 
data, LES, DNS or experimental definitively describing its stationary states. At the very least an 
assessment of the bifurcation diagram ( e/Sk ) 00 versus (£2/5) predicted by the linear theory would 
be useful. Also the equilibrium values ( V fe )oo would be useful for further developments regarding 
the dissipation equation’s dependence on rotation. This may remedy the under-prediction of the 
kinetic energy growth rates as a function of £2/ 5 for all the models. The present class of dissipation 
equations predicts an asymptotic state in which {V /s)oo a model dependent constant, which for 
all rotation rates has the same value that it does in the asymptotic shear. Had the stationary 
values of the anisotropy tensor and {V/e)oo an ^ (e/Sk ) oo been available application of the present 
methodology would have produced a set of modeled evolution equations whose fixed points matched 
to the fixed points of the rotating shear, independent of the deficiencies in the dissipation equation. 

The present rapid-pressure model is expected to distinguish itself in complex three-dimensional 
flows. For the planar flows for which test cases exist the model out-performs the nonlinear models 
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using realizability type constraints. It is only moderately better than the topologically generic 
form of the SSG model suitable for simple planar flows. It is unfortunate that there does not exist 
any suitable DNS or LES of flows in which the presence of a body force causes the larger scales 
of the motion to become quasi-two-dimensional. Such a test case would help further establish the 
utility of incorporating some of the more complex physics into the structure of the model as well 
as perhaps pointing out the potential deficiencies of models developed using two-dimensional mean 
flows. 

10. Summary and Conclusions 

A new representation for the rapid-pressure strain correlation with a minimum of ad hoc constants 
has been devised. The rapid-pressure model produces the proper behavior in five different limits: 

1) the geostrophic limit in which the eigenvalue of the Reynolds stress tensor, < UiUj >, along the 
axis of rotation vanishes, 

2) the realizable limit in which an arbitrary eigenvalue of the Reynolds stress vanishes, 

3) the joint-realizable limit in which an eigenvalue of < 96 >< U{Uj > — < 0u{ >< 6uj > vanishes, 

4) the isotropic limit in which the anisotropy tensor, b{j = 0, 

5) the asymptotic limit of a structural equilibrium in which DjDt b{j — 0. 

The model has the general form 

Xuu = Xgu + FXf ju + FXf° u 

where the X t jki are polynomials in the anisotropy tensor. X^- kl satisfies the five constraints: the 
limit states of geostrophy, realizability, joint-realizability, and the integral constraints of continuity 
and normalization. X^ kl satisfies the three constraints of joint-realizability, continuity and normal- 
ization. Both Xf jkl and FXf jkl are obtained analytically: they represent the simplest analytical 
expressions that are capable of satisfying all the mathematical constraints. X^ kl is obtained by 
requiring asymptotic consistency with the structural equilibrium state for homogeneous shear. Ex- 
perimental data is required to determine X^. The stationary values of the anisotropy tensor, 6^, 
and (V/e)oo and (e/Sk) oo are inserted into the modeled evolution equations to obtain the structural 
equilibrium component of the rapid-pressure, X^. This insures that the fixed points of the mod- 
eled equations match the experimentally determined fixed points. This process has been carried 
out for the homogeneous shear producing a model that should be accurate for three-dimensional 
mean flows in which the mean shear is the predominant production mechanism. 

Several points regarding the 2DMFI model merit mention: 

1) All the coefficients, A t , in the Xf- kl + FX[- kl portion of the model are obtained from first 
principles and are valid for all states of the mechanical turbulence. They take on their values at 
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the limit states only when the limit states are achieved unlike other models. 

2) The present form of the model is consistent with the stationary state of homogeneous shear. This 
insures that in the absence of any disturbing forces, the predicted flow will relax to the equilibrium 
state. 

3) The model is at present the only model which is consistent with predictions of the Taylor- 
Proudman theorem assuring that the modeled equations are frame-indifferent when the components 
of the Reynolds stress along the axis of rotation vanish. This frame invariance is the most notable 
feature of the model and is expected to be important for the computation of engineering and 
geophysical flows in which body forces play an important role. In meteorological flows the present 
2DMFI rapid-pressure correlation is consistent with the geostrophic limit atmosphere attains above 
the planetary boundary layer. As such the model will be relevant to the computation of mesoscale 
meteorological flows. 

It is expected that the satisfaction of the 2DMFI principle will be important in three-dimensional 
flows in which the largest scales of the motion are two-dimensionalized by body forces or kinematic 
constraints. These flows include: 1) turbulence in which a strong stable stratification suppresses 
the vertical component of the velocity field; 2) turbulence in the flows affected by magnetic fields 
3) turbulence influenced by rotation such as those in crystal growth processes or occurring in tor- 
nadoes, swirl combustors and turbines; 4) turbulence influenced by Coriolis forces in which the 
largest scales of the flow undergo a Taylor- Proudman reorganization (Ro < 2) as might occur in 
large scale geophysical flows; 5) turbulence near a free surface at which one of the components of 
the fluctuating velocity is suppressed and the mean flow is two-dimensional; 6) the environmen- 
tally important shallow water free shear flows as the shallow water jet associated with waste heat 
exchangers, near-shore pollution dispersal, tidal estuary flows, and mixing associated with thermal 
and salinity inflows in lakes and rivers. However, until suitable data bases, DNS or LES, of these 
complex flows with body forces become available the full potential of a rapid-pressure model built 
from first principles in three-dimensional flows can not be verified. The present model does however 
reproduce the experimental data at least as well as. the currently available models for a wide class 
of two dimensional flows. 
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Appendix 1: A synopsis of the final rapid-pressure correlation representation 

For convenience and clarity the final form of the model is summarized here. The general final form 
of the model is: 


IE; /2 q 2 = [C 3 - 2 IIC% + 3 IUC^Sij 

+ C 4 [b tp S PJ + bj p S p i - 2/3 < bS > 6 {j ] + C"' <bS> [&?• + 277/3*,,] 

+ Cs[bi p W p j + bj p W p i] + [C 6 <bS> +C"' < b 2 S >}b tJ 
+ C 7 [b 2 p S P] + b?S pi - 2/3 < b 2 S > 6 {j ] 

+ C s [b%W pj + b%W pi ] + C 9 [b ip W qp b 2 q] + b }P W qp b 2 qi }. 

For flows in rotating coordinate systems, Wij appearing in flij must be replaced by the total rotation 
tensor Wij + Cjik^k- The coefficients C t are given by a sum of the basic model coefficients, Bj, 
which come from first principles, and the calibration coefficients Aj. The basic model coefficients 
are related to the coefficients Ai in the fourth order tensor polynomial expression by 

B 3 = 2{A\ + A 2 ) = (2/27)[41 + 42/7 - 0.1F(221 + 42077)]/77 d 
= A 9 -+ A 10 = — (14/3)(1 + 3II)/II d + 0.6F/(1 + 377) 

B% = -1/3(A„ + A 12 + 2A 13 ) = (55 + 8477)/377j 
B 4 = A 3 + A 4 + 2 A 5 = 3/ 1 I d - 0.9F/(1 + 377) 

B' 4 — A\\ + Ai 2 + 4Ai3 = —9jII d 
Bs= A 3 -A 4 = — ( 1 / 30)( 10 + 21F)/(1 + 377) 

Be = 2A 9 + 4/lio = —1877/ 1 I d + 3F/(1 + 377) 

B r = A 6 + A 7 + 2 A s - 2A 9 - 2Aio = -9/II d - 1.8F/(1 + 377) 

B s = A 6 -A 7 = (1/5)(3F - 5)/(l + 377) 

B 9 = An — A\ 2 = —3/(1 + 377) 

The Aj are given in Appendix 3. The C, , which reflects the application of the additional constraint, 
asymptotic consistency with an equilibrium state, are obtained from the Bi by Ci = Bi + FAl 
They are given by 

C 3 = B 3 - 2F(10A| + 3 A c 9 + Af 0 )77/5 - F{A c n + A c 12 + \4A\ 3 ))III/b 

C% = B% + F(A% + i4f 0 ) 

Cf= B'z — 1/3F(Ajj + A\ 2 + 2A\ 3 ) 

C 4 = B 4 + F(— 3A| + II(A C U + A<j 2 + 4A\ 3 )) 

C? = B + F(Ajj + A c u + 4 A c 13 ) 

C 5 = B s + F(-7/3Al + (~A C U + 3 A c 12 + 4A\ 3 )II/3) 

C 6 = B 6 + F{2A% + 4A\ 0 ) 

C 7 — B 7 — 3F(A| + Ag + Aj 0 ) 

Cg = B 8 -l/3F(7A| + 3A|-Af 0 ) 

C 9 = B 9 + F(i4fj - A\ 2 ) 

where the Aj are expressed in terms of the seven free parameters A^{5,8— 13). Without further 
specifying the calibration coefficients, Aj, the above rapid-pressure model is general, suitable for the 
flows for which second-order modeling is suitable. Calibration to a particular archetypal flow, by 
matching the fixed points of the modeled evolution equations with experimentally determined fixed 
points, will result in a model suitable for diverse flows within that class of flows. For a flow in which 
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the dominant production mechanism is associated with the mean shear the hommogeneous shear 
is used to set the model coefficients. Asymptotic consistency with the homogeneous equilibrium 
shear and a linear model for the return coefficient produce the following values for the calibration 
coefficients: 

Ag = -0.29 - 0.06(Ag o - Ag) 

A c u = -3.6 + 5 A\ 0 - 2Aj 3 - 12.7Ag - 3.8Ag 
Ag 2 = -24.5 - 44.2Ai 0 - 2 Ag 3 + 29Ag - 8 Ag 

where Ag = 0.8, Ag = —1.0, Af 0 = 0.01, Ag 3 = 0. 

The homogeneous form of the continuity and normalization conditions have been used to produce 
relations between the free parameters. They are: 

A\ A 4 Ag - 2Ag// + III(A c n A A\ 2 + 2 A§ 3 ) = 0 
Ag A Ag A 5Ag - //(A^ + Ag 2 A 4 Af 3 ) = 0 
Ag A Ag A 5Ag + Ag A Af 0 = 0 

3Ag A 2Ag - 2Ag// + 4 A\ 3 III = 0 
3Ag A 4Ag - 2//(Af, A 2 Ag 3 ) = 0 
3Ag A 4Ag A 2Ag 0 = 0 

The fact that the Ag must satisfy these equations has been used to express six of the Ag{l —4,6, 7} 
in terms of the seven Ag{5,8 — 13} in the expressions for the C ,• = B± A FAg. 

Appendix 2 : The calibration coefficients in the rapid-pressure representation 

The model has been calibrated to match the fixed points of the homogeneous shear. The calibration 
process is described more detail in this appendix. The following equations describe the evolution 
of the anisotropy: 

DjDt bij = -2e,jfcp& p jftjt - 2e jkp b pi 9. k - [b ip Uj, p +b JP Ui, p -2/3 6 tJ < bS >] 

-2/35, j A 2 bij < bS > A EG,- /<? 2 - (Cl - 2 )b t] e/q 2 A C 2 [b% A 2II/36 ij ]e/q 2 

For the planar flow case with axial mean flow and mean shear, t/ 1 , 2 , the algebraic equations for the 
fixed points, setting D/Dtbij = 0 , become 

(b n A l/3)(2 A 4 b u Sk/e) - 4b 12 Sk/e A 4 (Skfe) Fin /(4 kS) - Ci 6 „ - 2/3 
A C 2 (b u bu A & 12&12 A 2/7/3) - 8(Cl/ S)(Ske)b u = 0 

612(2 A 4b\ 2 Sk j e) — 2(622 A l/3)5fc/f A 4 (Sk/ £) FU 2 /(4&5) — C\b\2 
4-C 2 b\ 2 (bn A 622 ) A 4(Q / S)(S k / e)(bi\ — 622 ) = 0 

(622 A l/3)(2 A 4b\ 2 Skj s) A 4(Sk/s') O 22 /(4fc5) — C\b 22 — 2/3 
+£ 2(622622 + 612612 A 2/7/3) A 8({l/S)(Sk/e)b u = 0 

(633 A l/3)(2 A 4b\ 2 Skje) A 4 (Sk/e) ng 3 /(4 5fc) - £^33 - 2/3 A £ 2(633633 + 2///3) = 0 

The last equation for 633 is not linearly independent as bjj = 0 . The following fixed point values are 
taken to be representative: b^ = 0.203, 6^2 = -0.156, 6 ^ = -0.143, 633 = -0.06, (5A:/e ) 00 = 5.54, 
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(P/e ) oo = 1-73. are inserted into the algebraic equations above. The values of the invariants 
corresponding to these values of the anisotropy tensor are: 7 /qo = —0.0604, 1 1 loo = 0.0038, Foo = 
0.5613. The solution of the algebraic equations describing the stationary state of the homogeneous 
shear produces the following values for the calibration coefficients, 

A% = -0.3388 - 0.06(j4g + A\ 0 ) + 0.015C? 0 - 0.008(72° 

A c u = -14.35 - 2A c 13 - 12.68 A% - 3.80A| - 5.072^ o + 3.157(7 ,°° + 0.5898<7 2 °° 

A c u = -16.87 - 2 A\ 3 + 28.7Ag - 8.05^1g - 44.81A* 0 - 12.17 <7f - 0.291(72°. 

The model has been left as general as possible - allowing for any return term of the canonical form, 
-C\sbij + C' 2 [6 + 21 1 /3<$,j], were the (7, are not necessarily constants but have achieved their 
asymptotic values <7, = Cf°. Choosing a linear Rotta type return term setting the nonlinear return 
coefficient to zero, C 2 = 0, produces 

A% = -0.29 - 0.06(.4j O - A c a ) 

A c u = -3.6 + 5j4j 0 - 2A C 13 - 12.7Ag - 3.8 A% 

A c u = -24.5 - AA.2A C W - 2 A\ 3 + 29Ag - 8 A c 9 

where the following values of the free parameters have been chosen: A\ = 0.8, Ag = — 1.0, A c l0 = 
0.01, A\ 3 = 0. These indicated calibration coefficients are appropriate for the class of flows in which 
the mean shear is the predominant production mechanism. 

The component form of the rapid-pressure used in the equations for the stationary state is 

Fin /2g 2 5* = 
n r 2 2 /2g 2 s* = 
m 2 /2 q 2 S* = 


n r 33 /vs* = 


C4612/3 + <T(&n&ii + 612612 + 2 /// 3 ) 6 i 2 + C 5 b 12 (-W*/S*) 

+[(76 + C'/'{b\\ + 622)1611612 + <77(611 + 622)612/3 

+< 7 8 ( 6 u + b 22 )b l 2 {-W*/S*) + C 9 (W*/S*)b u (b 12 b l2 - 611622) 

C4612/3 + C !/' (622622 Hh 612612 + 2 /// 3 ) 6 i 2 + C$bi 2 (\V / 5 ) 

+[< 7 $ + < 74 ,, ( 6 ii + 622)1622612 + <77(611 + 622)612/3 

+< 7 8 ( 6 h + 622)612(^/5*) + C 9 (-W*/S*)b l 2 (b 12 b 12 - 611632) 

l/ 2 [( 7 3 - 2//C3 + 3///(73'] + l/ 2 ( 7 4 ( 6 n + 622) + C'"(6n + 622)612612 
+<75(611 - b 22 )(W*/ 2 S*) + [< 7 6 + CT(6ii + 622)1612612 
+ 1/2(77(6 h 6 ii + 2612612 + 622622) 

+< 7 8 ( 6 ii 6 h — b 22 b 22 )(W* / 2 S*) + C 9 (bnb 22 — 6i 2 6i2)(6n — b 22 )(W* / 2 S*) 

— 2/3C4612 + C'/^b^b^ + 277 / 3 ) 6 i 2 + (<76 + < 74 "( 6 n + 6 2 2 )) 6 i 26 3 3 
—2/3(77612(611 + 6 22 )- 
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Appendix 3: The representation of the integral of the two-point covariance 

The general form of the fourth-order tensor polynomial used to model the volume integral of the 
derivative of the two-point covariance, X x jku appearing in the Reynolds stress equation is 


Xijki/ < u p u v >= A\f*ij8ki + A2(6ik8ji + SuSjk) 

+ A 3 6ijbki + A 4 bij8ki + A&(bik6ji 4- + t>ikbji + Subjk) 

+ A^6{jb^ + ArtifjSki + A&(tif k 6ji + b^Sjk + 8%kb*i + bubj k ) 

+ A 9 bijbkt + Aio(bikbji + bubjk ) 

+ Anbijb^ + Awtifjbki + A\$(b* k bji + b^bjk + + bub^ k ) 

+ 

The coefficients with all free parameters set to zero, as derived by Ristorcelli (1991) or in Ristorcelli 
and Lumley (1987,1991), are 

Aj = (111/7+ 73)/27// d - F(420// + 239)/135//<i 
A 2 = -(69/7 + 32)/27// d + F(420/7 + 257)/27077 d 
A 3 = (377 + 4)/3// d - F(ll/10)/(1 + 377) 

A 4 = (1577 + 11 )/ZII d - F(4/10)/(l + 377) 

A 5 = -3(1 + 377)/377 d + F(3/10)/(l + 377) 

A 6 = -(10277 + 61)/377 d 

At = -2(3377 + 20)/3// d - F(6/10)/(l + 377) 

+ 8 = (4277 + 23)/377 d 

Ag = -(5777 + 28)/377(i - F(3/10)/(l + 377) 

Aio = (1577 + 14)/377rf + F(9/10)/(l + 377) 

An = -(10277 + 61)///* 

^12 = -2(33/7 + 20)///* 

A 13 = (42/7 + 23)///* 

Al4 = 0 

Au= 0 

where 77* = (l+377)(7+1277), F = 1 +27777+977, where 77 = —1/26^6^ and 777 = 1/3 b ip b pj bji. 
The rapid-pressure correlation integral, X p kj , appearing in the heat-flux equations, is used to 
derive the A{ through the joint-realizability constraint. The general form of the third-order tensor 
polynomial used to model X p kj is included for completeness: 

Xpkj — 71 1 0u p 6k j + D 2 [< Silk ^ ^ Ouj ^ 

+ D 3 < 6u p > + 7?4[< flu* > b p j+ < 0Uj > ftp*;] 

+ D 5 < 6u p > b 2 kj + £> 6 [< 9u k > bpj+ < 9uj > b 2 pk ] 

+ [DjbqpSkj + D$(b q k6 p j + b q j 6 pk )] *C 9u q > 

+ \Do>bqpbkj + Dio(b q kbpj + b q jb p k )] + 0u q > 

+ [Diibq P b 2 j + Di 2 (bqkb p j + b q jb pk )] < $u q > 

+ [D\zb 2 p 6kj + Di4(b 2 k 6 p j + b 2 j6 p k)] < Qu q > 

+ [Disb 2 qp b k] + Die(b 2 b PJ + 62.6 pfc) ] < Q Uq > 

+ [Z? 17 62 p6 2 . + Di 8 (b 2 k b 2 pj + b 2 jb 2 k )] < 0u q > 


34 


Where the D, are given in Ristorcelli and Lumley (1987,1991) as 


£>i = -(312 II 2 + 1497/ - 21)/5 II d - F( l/5)/(l + 3/7) 

D 2 = (4877 2 + 77 - U)/bII d + F(3/10)/(l + 377) 

D 3 = — (32477 2 + 22277 + 17)/77,i - F(3/10)/(l + 377) 

D 4 = -3/(7+ 1277) + F(9/10)/(l + 377) 

D 5 = -(10277 + 61)/77 d 

D 6 = (4277 + 23)/77 d 

D 7 = -2(377 + 4)/577 d - F( 3/5)/(l + 377) 

7?8 = 27(277 + l)/577d 
D 9 = (4277 + 23)/ 1 Id 
Dio = (4277 + 23)/ 1 I d 
D 13 = -8(3977 + 22)/5IId 
D\4 — 2(2477 + 17)/ 577^ 

£>15 - -27/(1 + 377) 

and Du = D12 = Di6 = D17 = Dis = 0. The first term of each of the coefficients 4,-, Di satisfies 
the constraints of continuity, normafization, realizability, joint-realizability and 2DMFI. The terms 
proportional to F are the terms that contribute to the rapid-pressure when the turbulence is away 
from both the realizable and the 2DMFI limits. 

Appendix 4: The Cayley-Hamilton theorem generalization 

Reference has been made to a generalized Cayley-Hamilton theorem, Rivlin (1955), relating different 
powers of products of matrices. 

ABC +ACB + BCA + BAC + CAB + CBA = 

A{< BC > - < B >< C >) + B{< CA > - < C >< A >) + C(< AB > - < A >< B >) 
+{BC + CB) < A > +(CA + AC) < B > +{AB + BA) <C> 

+ 1[< A >< B >< C > - < A >< BC > - < B >< AC > - < C >< AB > 

+ < ABC > + < CBA >] 

The Cayley-Hamilton generalization is easily derivable from the Cayley-Hamilton theorem applied 
to sums and differences of matrices A, B,C. The theorem is useful in eliminating redundant tensor 
bases in tensor representation theorems. Here <> is used to indicate the trace. The Cayley- 
Hamilton theorem for the anisotropy tensor is b 3 = 1/3 < 6 3 > 1 + 1/2 < 6 2 > b. The theorem can 
also be used to express 

bSb = — [b 2 S + Sb 2 ]+ < bS > b + 1/2 < b 2 > S+ < b 2 S > 1 
bSb 2 + b 2 Sb = -1/3 < b 3 > S+ < b 2 S > b+ < bS > b 2 

Appendix 4: The rapid-pressure' models 

The detailed form of the pressure-strain models referred to in this paper are as follows: 

The Launder, Reece & Rodi Model 

4 2 

Uij = —2 Ci£b {j + -KSij + C 2 K(b tk S jk + b jk S ik - -b k ,s ki 6 t: ) + C 3 K{b tk W }k + b jk W ik ) 

5 4 
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where 

Ci = 1.5, C 2 = 1.75, C 3 = 1.31 


The Shih & Lumley Model 

Hij = -febij + - A' ,5V, + 12 a 5 K(bikSjk + bj k Si k --b k iS k [6ij) + -(2 - 7as)K(bi k Wj k + bj k Wi k ) 

4 4 

"I - ^ K {bubhnS jjri 4* bjibimSim 2b lk S fc/ bij 3bkiSkibij) 4" ~^F(bnbi m Wj m 4“ bjibi 7n Wj rn ') 

where 

(3 = 2 4- j exp(-7.77/v/^){72/^Ae^ + 80.1 ln[ 1 4- 62.4(-/7 + 2.3 III)}} 

F = 1 + 97/4- 27777, 77 = OI ® \b tJ b ]k b kt 


„ 4 A' 2 1 / 4 i \ 

flet= 9^> Q5 = To( 1 + 5 F V- 


The Fu } Launder & Tselepidakis Model 

1 4 2 

Ilij = (3\eb,j + /3 2 e(b ik b k j - -b k ib k ib XJ )+-K S l} 4- l.2K(b lk S jk + bj k Si k - - b k iS k [6ij ) 

26 4 

+—K(b ik Wj k 4- bj k W ik ) + -K(b ik b k iSji+bj k b k iSu - 2b ik S k ib\j - 3b k iS k ib tJ ) 

+ ^K(b ik b kl W J , + b ]k b kl W tl ) - ^-K[8II(b tk W jk + b jk W ik )+12(b ik b kl W lm b mj + b jk b kl W lm b mi )] 
o o 


where 


ft = 12077F 1 / 2 4- 2 F x < 2 - 2, ft = 144//F 1/2 


The Speziale, Sarkar & Gatski Model 

n tJ = -(2Ci£ + c;v)bij + C 2 e (b ik b kj - i, bubuti ,•) 4- (C 3 - C 3 1 1 If )K Si j 
+C 4 K (b ik S jk 4- bjkSik - -bkiSkitiij'j 4- C 5 K(bi k Wj k 4- bj k W tk ) 

where 

Ci = 1.7, C* - 1.80, C 2 = 4.2 

C 3 = 1 C 3 * = 1.30, C 4 = 1.25 
5 

Cs — 0.40, lib = bijbij 
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Figure 3. Off-diagonal component of the anisotropy tensor as a function of the stability 
parameter for the homogeneous shear with curvature: 2DMFI , SSG , SL 
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Figure 4. Bifurcation diagram for the homogeneous shear with curvature for the four different 
models. 2DMFI , SL , and FLT , SSG . 
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Figure 5. Time evolution for the anisotropy tensor in plane strain for eo/Sko = 2.0. The 

predictions of the 2DMFI model , the SL model , and the direct numerical 

simulations of Lee and Reynolds (1985). 
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Figure 6. Time evolution for the anisotropy tensor in the axisymmetric contraction for 

eo/Sko = 0.179. The predictions of the 2DMFI model , the SL model , and the 

direct numerical simulations of Lee and Reynolds (1985). 
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Figure 7. Time evolution for the anisotropy tensor in the axisymmetric expansion for 

eo/Sko = 2.45. The predictions of the 2DMFI model , the SL model , and 

the direct numerical simulations of Lee and Reynolds (1985). 
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